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Abstract. In the present paper we generalize the original work of C.W. Misner 1 1 ] about the quantum dynamics of the Bianchi 
type IX geometry near the cosmological singularity. We extend the analysis to the generic inhomogeneous universe by solving 
the super-momentum constraint and outlining the dynamical decoupling of spatial points. Firstly, we discuss the classical 
evolution of the model in terms of the Hamilton- Jacobi approach as applied to the super-momentum and super-Hamiltonian 
constraints; then we quantize it in the approximation of a square potential well after an ADM reduction of the dynamics with 
respect to the super-momentum constraint only. Such a reduction relies on a suitable form for the generic three-metric tensor 
which allows the use of its three functions as the new spatial coordinates. 

We get a functional representation of the quantum dynamics which is equivalent to the Misner-like one when extended point by 
point, since the Hilbert space factorizes into o° 3 independent components due to the parametric role that the three-coordinates 
assume in the asymptotic potential term. 

Finally, we discuss the conditions for having a semiclassical behavior of the dynamics and we recognize that this already 
corresponds to having mean occupation numbers of order ff(W 2 ). 
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1. INTRODUCTION 

The evolutionary scheme of the universe is well described by the Standard Cosmological Model (SCM) from the 
light-elements nucleosynthesis (~ 10~ 3 — 10~ 2 sec.) up to the present stage of the observations; however, the very 
early phases of the cosmological dynamics (i.e. above 10 14 — 10 15 GeV) requires more general schemes for being 
completely understood. 

In fact, either the backward instability of the isotropic universe with respect to tensor perturbations either 
the impossibility to preserve symmetry in quantum gravity, strongly support the idea that the early universe must be 
represented by a generic inhomogeneous model, as discussed in jj,|5|]. 

As well known, the point-like dynamics of such a generic picture closely resembles the one characterizing the 
Bianchi type VIII and IX geometries, i.e. the so-called Mixmaster |6]. First investigations about the classical and 
quantum behavior of the Mixmaster were provided by C.W. Misner in |6] and 1 1], respectively; his analysis relies on 
a Hamiltonian formulation of the dynamics and on a suitable approximation of the spatial curvature which here plays 
the role of a potential term for the generalized coordinates (universe volume and anisotropics). 

The present work aims to generalize the main features of the Misner classical and quantum approaches toward the 
generic inhomogeneous universe, by virtue of the dynamical decoupling characterizing spatial points near the Big 
Bang. 

The classical inhomogeneous dynamics is investigated in Section [2] via a Hamilton- Jacobi approach to the funda- 
mental constraints of the theory; the solution outlines the right number of four physically arbitrary functions of the 
spatial coordinates in the configuration space as required by the generality. 

In Section[3] before proceeding to the quantization, we perform an Arnowitt-Deser-Misner (ADM) reduction of the 
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super-momentum constraints; in this way, we remove any dynamical role of the spatial gradients from the quantum 
behavior. Indeed, asymptotically to the singularity, the three coordinates enter the potential term as parameters only 
and then we can apply the long-wavelength approximation. Such a framework has the physical meaning of dealing 
with inhomogeneous scales which are super-horizon sized; thus, when below we will speak about spatial points, we 
will refer to causal regions of the universe which asymptotically are dynamically decoupled. This fact provides a sort of 
causal structure to our quantization procedure in the sense that no-interference phenomena take place on super-horizon 
scales. 

In Section |4] we give a quantum representation of the dynamics in terms of a functional approach which outlines 
the factorization of the Hilbert space into °° 3 point-like components. The quantum dynamics relies on an adiabatic 
approximation ensured by the potential term behavior, which reduces to an infinite well; according to C.W. Misner 
we model the potential as an infinite square box having the same measure of the original triangular picture. The 
volume-dependence of the wave functional acquires increasing amplitude and frequency of oscillations as the Big 
Bang is approached and the occupation number grows, respectively. 

Finally, in Section [5] we discuss the conditions under which the WKB semiclassical limit is recognized. We show 
that this corresponds to having occupation numbers sufficiently greater than unity and therefore semiclassical packets 
of eigenstates would correspond to average mean occupation number of order ^(lO 2 ). 



2. GENERIC COSMOLOGICAL MODEL 

When considering a cosmological solution containing a number of space functions such that a generic inhomogeneous 
Cauchy problem is satisfied on a non-singular hypersurface, we refer to it as a generic inhomogeneous model 
In the Arnowitt-Deser-Misner (ADM) formalism, the corresponding line element reads as 



ds 2 = -N{t,x) 2 dt 2 + Yijidx 1 +N i dt)(dx j +N j dt) , (1) 

jet ja 



where N is the lapse-function, N' the shift- vectors and q a (t.x') three scalar functions (the index a is summed over 
1,2,3); the vectors If have components which are generic functions of the spatial coordinates only, available for the 
Cauchy data. The general case, in which /" are time-dependent vectors is addressed in a following Section to simplify 
the variational principle. It is convenient to introduce also the reciprocal vectors l' a , such that lfl' b — 8g and I" l' a = 5j. 
The dynamics of this system is summarized by the action 



[ d i xdt{p a d,q a -NH-N i H i \ (3) 
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H i = -2 X d j 'Eliirp m -XP } diqj, (5) 
m 

being p" the conjugate momenta to the variables q a , y = det^-, the tridimensional Ricci scalar and % the Einstein 
constant. By variating the action Q with respect to the N and N', we get the super-Hamiltonian and super-momentum 
constraints H = and //, = 0, respectively. 

Let us introduce the Misner-like variables via the transformation 

qi = 2a + 2/3+ + 2 V3j3_ (6a) 
q 2 =2a + 2/3+-2V3/3_ (6b) 
g 3 = 2a -4/3+; (6c) 



in terms of this set of configurational coordinates the Hamiltonian constraints rewrite as 

V 
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A detailed analysis of the potential term V leads to [2] 



V =t- {af(jc'>- 8 ^ +al(x i )e 4 ^ + ^+al(x i )e 4 ^-^+W(x i , a, j8±, dja, djP ± )} , 



where a, (i — 1,2,3) refer to the space quantities 



a i (jr / )==r'-rotl i , 



and we regard the operations • and rot as taken in Euclidean sense. 

To outline the relative behavior of the two terms in the potential as the singularity is approached for a 
consider the quantities 
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Taking into account these definitions, the potential V rewrites as 
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Near the cosmological singularity a — > — °° and D — > 0, so that the term W becomes negligible. Indeed this conclusion 
is supported by the behavior of the spatial gradients, which does not destroy the feature above outlined (see 0,01). 
Through the canonical replacements 
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the classical evolution is summarized by the Hamilton-Jacobi system 
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Since sufficiently close to the cosmological singularity the potential term becomes negligible, then the solution of 
illai reads as 



k 2 + + k 2 _ a + £ + j3+ + /t_/3_, (18) 

where k± — k± (x l ) are arbitrary functions of the coordinates and the minus sign in front of the square root has been 
taken because we are considering an expanding Universe. 

According to the Jacobi prescription, the functional derivatives of the above action Jl 8I > with respect to k± have to 
be set equal to stationary quantities c± (x l ) and therefore we get the following expressions for j3± in terms of a 

P± = 7t ± a, (19) 

where 

k± 

7C± = (20a) 

lkl+kl 
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Substituting the solution Jl 8i with dl9l in the Hamilton-Jacobi equation J17bl corresponding to the super-momentum 
and taking into account the relations j20l . the last sum cancels out leaving the equations 



-2 ++ di(k + +k-) + dik + + V3dik-+dj \l J 3 lf (k+-2V3k-) -2V3l J 2 lfk- 
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= 0, (21) 



which are constraints on the spatial functions only. The above mentioned functions c±(x l ) have been set equal to zero 
because their presence would simply correspond to a rescaling of the vectors lf(x J ). Thus, our solution contains ten 
arbitrary functions of the spatial coordinates, namely the nine components of the vectors If and one of the two functions 
7t±. These ten free functions have to satisfy the three constraints Hit : the choice of the coordinate frame eliminates 
the arbitrariness of three more degrees of freedom. Therefore our solution is characterized by four physically arbitrary 
functions of the spatial coordinates and, in this sense, it is a generic one. 

Above we neglected the role of the potential because it influences the point-Universe evolution only via the bounces 
producing the establishment of a new free motion (for a detailed discussion about the chaotic properties of the random 
behavior that the point-universe performs in the potential, see l9l [l(il.llTlo . This effect of the potential is summarized 
by the reflection law 

sin 0/ - sin 0,- = - sin(0/ + 0,-) , (22) 

where 0,- and Of denote the angles of incidence and deflection, respectively, of the point-Universe for a bounce on one 
of the three equivalent walls of the triangular potential V(a,j3+, j3_), see Figure^ 

The results of this Section show how the generic cosmological solution toward the Big-Bang is isomorphic, point 
by point in space, to the one of the Bianchi types VIII and IX models lfl2ll because the spatial coordinates are involved 
in the problem only as parameters. 



3. ADM REDUCTION OF THE DYNAMICAL PROBLEM 

In the last Section we have discussed a generic cosmological model in which the vectors If appearing in (0 were 
functions of the coordinates x ; only. It is interesting to discuss a more generic framework, that is to allow the more 
general dependence If — lf(t,x } ). 

To this purpose, let us rewrite such vectors in the form fl3ll 

lf = O a b d i y b , (23) 

where O a b = Of / (x k ) is a SO(3) matrix and y h — y b (t,x k ) are three scalar functions. Considering the metric tensor (|2ji 
rewritten as 

Yij = e^8 ad O a b O d c d i y b d j y c , (24) 
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FIGURE 1. Equipotential lines showing the potential domain Th , dynamically closed in the corners: this visualization is in the 
j6+, j6_ plane of the Misner variables. Asymptotically, the walls flatten to a triangular well, as in Equation 1301 . 



with the indexes (a,b,c,d,i,j) running as 1,2,3, the action (|3}becomes 



where, apart some details of the potential coordinate dependence, H has the same expression as Q, while the 
super-momentum explicitly reads 



Yl c i being the conjugate momenta of the variables y d . 

Summarizing, we have ten independent functions characterizing the dynamical system: the 3 degrees of freedom of 

the coordinate choice y", 3 scale factors q", 3 components of the shift vector N' and the lapse function N. 

The variation of the action fl25i with respect to N and N' leads to the usual constraints H = and //, = 0, respectively. 

Choosing y" as spatial coordinates, i.e. y" — y"(t,x) and rj = t, we can diagonalize and solve the constraint d26l as 

follows 
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Given |7| as the Jacobian of the change of coordinates, the relations 



< 




(28) 



, dx a dx a dy b 



allow to write the (reduced) action as 
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being 



«*>- C -t"^ (31) 

if x < . 



The dynamics of the universe is decoupled for each space point and the point-universe moves in the domain r# as 
outlined by the triangular equipotential lines of Figure^ which is dynamically closed also in the corners. 

The vanishing behavior of the potential inside the domain of definition implies that, near the singularity, dp a /dr\ = 
0. Therefore the term p a (0~ l ) d^O^. in d29i reads as an exact time differential term, and unaffects the asymptotic 
dynamics whose variational principle takes the form l i 1 311 



RED — 



[ n d7]d 3 y{p a d n q a -NH), (32) 



and the dynamics is described by H only. 



4. GENERIC QUANTUM COSMOLOGY 

Let us now analyze the quantum dynamics corresponding to a generic cosmological model within a canonical 
framework, by replacing the canonical variables with the corresponding operators and implementing the Hamiltonian 
constraints on the state functional describing the system, *P = *P(oc, j8+,j3_). Therefore we adopt the representation 
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-ih- 



'8a' Pp± ' "'S/3±' 
and implementing H to an operator the quantum dynamics is described by the Wheeler-deWitt equation 
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(33) 
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As discussed at the end of the preceding Section, since the asymptotically the potential term remains isomorphic to 
the Bianchi IX one, point by point in space, we can find a solution in the form 



V/ = £r„(a)<p„(a,j3+,j3_) : 



(35) 



where the coefficients F„ 



r„ =exp / d 3 xlnC n{xi) (a) = Y\C n{xl) (a 



keep trace of the inhomogeneities of the system and the index n has to be regarded as a space function n(x l ). 
Thus, equation ( 1341 is reduced to the ADM eigenvalue problem 
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According to Q| we approximate the triangular infinite walls of the potential by a box having the same measure and 
then we find the eigenvalues E n 



M \ ' \n\ 
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(38) 



where n 2 = n 2 ,+n 2 _, being n + ,n_ 6 ,JV the two independent quantum numbers corresponding to the variables j3+,/3_ 
respectively. 

Substituting the expression for F„ in equation j34l we get the differential equation 



Y,(d 2 C n )(p„ +£C„(<?2 (pn ) +2Y,(d a C„)(d a <p„) +Y,E 2 C„<p n = 0, 



(39) 




which, in the limit of the Misner adiabatic approximation of neglecting d a <p„ (i.e. <p ~ 0(j3 + ,/3_)), simplifies to 

7 d 2 C n k 2 

h ^ + ^ C "=° (40) 



where 
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%=[^\ \n\ 2 h 2 . (41) 



The above equation is solved by C„(a) in the form 

C„(a) = CiVasin [ i^/p^lna ) + C2\/acos ( ^-yfp^lna ) , (42) 



where y/p^ = ^Jk 2 — 1. Figure |2] shows the behavior of C„(a) for various values of the parameter k n . Such wave 
function at fixed x' behaves like an oscillating profile whose frequency increases with occupation number n and 
approaching the Big Bang, while the amplitude depends on the a variable only. 



5. SEMICLASSICAL LIMIT OF THE QUANTUM DYNAMICS 

It is now interesting to consider the evolution of the system in the semiclassical limit, that is to investigate the dynamical 
behavior in the vicinity of the Big Bang within the WKB approximation of the dynamics (for a discussion of the 
semiclassical limit in the homogeneous case, see Q). To this purpose, let us search for a solution of the wave 
equation J40i in the form 

C(a)-Vp e ,ff / fi , (43) 



which reduces such equation to the system of the real and complex parts as 



respectively. In the semiclassical limit, neglecting the term in square brackets in Equation J44a> . we obtain the solution 



a(a) = klna, 
p(a) = const, x a. 

This behavior remains valid as far as the term in h 2 is negligible and therefore we get the constraint 



(45a) 
(45b) 




(46) 



i.e. n 1. As soon as this constraint is fulfilled, the solution remains semiclassical close to the Big Bang and we 
see that, in the sense of the WKB approximation, this would imply packets of eigenfunctions whose mean quantum 
number <n> can be of order ^"(lO 2 ) only. 



By our analysis, we gave an inhomogeneous extension of the classical and quantum Mixmaster dynamics. Following 
the original point of view of C.W. Misner, we achieved a quantum description of the model by approximating the 
potential term as an infinite square box. 

The main issue of our investigation has shown how the removal of the spatial gradients dynamics is possible in the 
asymptotic regime. This result comes out from an ADM reduction of the super-momentum constraints and from the 
parametric role that the three-coordinates assume in the potential term. 

On the quantum level, this achievement turned out into the factorized structure of the Hilbert space; from a physical 
point of view, this implied to deal with independent quantum behaviors within each cosmological horizon. 

Finally, we discussed the conditions allowing a WKB semiclassical approach to the dynamics and we fixed as 
condition to have occupation numbers much greater than unity. 
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6. CONCLUDING REMARKS 
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